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CONFERENCE  PAPER 


Recent  research  has  highlighted  the  need  for  physically  consistent  radiation  pressure  and  Bidirectional  Reflectance 
Distribution  Function  (BRDF)  models.  This  paper  seeks  to  evaluate  the  impact  of  B RDF -consistent  radiation  pres¬ 
sure  models  compared  to  changes  in  the  other  BRDF  parameters.  The  differences  in  orbital  position  arising  because 
of  changes  in  the  shape,  attitude,  angular  rates,  BRDF  parameters,  and  radiation  pressure  model  are  plotted  as  a 
function  of  time  for  simulated  rocket  bodies  and  debris  at  geosynchronous  orbit  (GEO).  The  initial  position  and 
velocity  of  the  space  object  is  kept  fixed,  and  the  orbital  position  difference  between  a  baseline  orbit  and  the  per¬ 
turbed  orbit  are  plotted  as  a  function  of  time.  This  is  similar  to  how  the  effects  of  perturbations  have  been  visualized 
in  the  past  in  commonly  used  astrodynamics  references. 


1.  INTRODUCTION 

There  are  many  sources  of  perturbations  when  propagating  an  orbit  forward  in  time  including  non -spherical 
gravitational  effects,  third-body  effects,  and  radiation  pressure.  Solar  radiation  pressure  (SRP)  depends  on  both  the 
space  object  (SO)’s  attitude  and  the  surface  properties  as  expressed  by  the  BRDF.  Wetterer  et  al.^  first  demonstrated 
how  the  SRP  acceleration  should  be  physically  consistent  with  the  associated  BRDFs  to  achieve  high-precision  orbit 
propagation. 

The  work  presented  here  seeks  to  compare  the  magnitude  of  the  perturbations  caused  by  using  a  simplified  SRP 
acceleration  model  as  compared  to  the  physically-consistent  SRP  model  with  those  resulting  from  changes  in  other 
parameters  associated  with  the  SRP.  The  initial  position  and  velocity  of  the  space  object  is  kept  fixed,  and  the  or¬ 
bital  position  difference  between  a  baseline  orbit  and  the  perturbed  orbit  are  plotted  as  a  function  of  time.  This  is 
similar  to  how  the  effects  of  perturbations  have  been  visualized  in  the  past  in  commonly  used  astrodynamics  refer¬ 
ences^.  First,  the  shape  model,  BRDF  model,  and  the  SRP  model  are  described.  Finally,  the  simulation  setup  and 
results  are  presented. 
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Figure  1.  Shape  Models 
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2.  SHAPE  MODEL  DESCRIPTIONS 


Each  shape  model  is  composed  of  a  number  of  flat  facets  with  each  facet  having  a  surface  area,  position,  and  nor¬ 
mal  direction.  Six  shape  models  are  considered  in  this  paper.  Figure  1  shows  a  rendering  of  each  model  and  Table  1 
displays  relevant  information  about  each. 

Table  1.  Shape  Model  Descriptions 


name _ mass  (kg) 


simplerocketbody 

1500 

lUSrocketbody 

300 

SLllrocketbody 

2400 

simpledebris(  1 ) 

1 

simpledebris(2) 

5 

crinckleddebris 

5 

_ description _ 

9  X  3  m  (cylinder)  -  rough  mass  and  dimensions  of  Cen¬ 
taur  upper  stage 

1.5  X  2.0  m  (truncated  cones)  -  rough  mass  and  dimen¬ 
sions  of  Inertial  Upper  Stage  (lUS)’s  stage  II  and  spacecraft 
interface 

4.9  X  3.7  m  (truncated  cones)  -  rough  mass  and  dimen¬ 
sions  of  Blok  DM  SL-12  upper  stage 
2  X  1  X  0.1  m  (rectangular  prism) 

2  X  1  X  0.1  m  (rectangular  prism) 

2  X  1  X  0.1  m  (rectangular  prism  with  49  facet  on  ±  z- 
_ faces  with  U  random  spread  of  normal  vectors) _ 


3.  BIDIRECTIONAL  REFLECTANCE  DISTRIBUTION  FUNCTION  MODEL 

Each  surface  on  the  shape  model  may  represent  a  different  material,  and  each  material  that  makes  up  the  model 
could  reflect  light  differently.  Thus,  a  model  might  have  many  parameters  that  specify  its  surface  properties.  The 
function  that  defines  how  light  is  reflected  from  an  opaque  surface  with  a  given  surface  normal  direction  (N), 

illumination  direction  ( L  with  angles  ft  and  from  N  ),  and  observer  direction  ( V  with  angles  ft  and  from  N )  as 
shown  in  Fig.  2  is  called  the  bidirectional  reflectance  distribution  function  (BRDF). 


The  BRDF  is  given  by 


Figure  2.  The  Geometry  of  Reflection 


dLXe„(t>/,X) 
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(1) 


where  dLr  is  the  reflected  radiance  in  Wm'^sr^  and  dEi  is  the  irradiance  in  Wm'^.  The  bisector  vector  between  the 
illumination  source  and  the  observer  is 


with  angles  a  and  ^  from  N  and  is  used  in  many  analytic  BRDF  models. 

There  are  many  different  reflectance  models  that  could  be  used,  but  all  can  be  expressed  in  a  common 
nomenclature  with  the  general  BRDF  calculated  by 

fr={dR^+sR^)  (3) 


which  depends  on  the  diffuse  bidirectional  reflectance  {Rd)  and  the  specular  bidirectional  reflectance  {Rs)  and  the 
fraction  of  each  to  the  total  (J  and  5*  respectively  where  d  s  =  \).  These  bidirectional  reflectances  are  calculated 
differently  for  the  various  models.  In  this  paper  we  will  use  the  Ashikhmin-Shirley  BRDF^,  also  known  as  the 
Anisotropic  Phong  BRDF,  where  the  diffuse  and  specular  bidirectional  reflectances  are  calculated  using 
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where  (4)  is  a  non-Lambertian  diffuse  BRDF,  and  the  Fresnel  reflectance  (F)  in  (5)  is  given  by  Schlick’s 
approximation"^ 


F  = 


(6) 


In  addition  to  J,  p  and  Fo,  the  Ashikhmin-Shirley  BRDF  has  two  exponential  factors  {Uu,  Uv)  that  define  the  ani¬ 
sotropic  reflectance  properties  of  each  surface.  Without  loss  of  functionality,  the  diffuse  reflectance  and  the  specu¬ 
lar  reflectance  at  normal  incidence  can  be  set  equal  to  each  other  (p  =  Fo)  and  the  difference  between  the  diffuse  and 
specular  reflectances  displayed  in  the  diffuse  fraction  parameter,  d.  Additionally,  for  the  sake  of  simplicity,  in  this 
paper  the  two  exponential  factors  are  set  equal  to  each  other  as  well  (riu  =  nv  =  n).  Thus,  there  are  three  unique  sur¬ 
face  parameters  per  surface  {n,  p,  d).  These  surface  parameters  have  constraints  (^>0,  0<yO<l,0<(i<l).  A  final 
assumption  is  that  every  facet  on  the  shape  model  is  represented  by  the  same  BRDF  parameters.  In  the  simulations 
of  the  next  section,  the  reflectance  parameter  is  set  to  p  =  Fo  =  0.1,  0.5,  and  0.9  representing  low  reflectivity,  mid 
reflectivity,  and  high  reflectivity,  the  diffuse  fraction  parameter  is  set  to  J  =  0.2  and  0.8  representing  a  highly  specu¬ 
lar  and  a  highly  diffuse  surface,  and  the  exponential  factor  parameter  is  set  to  =  ^  =  50  and  500  representing 

a  “dull”  and  “sharp”  specular  peak. 

The  curved  cylindrical  surface  of  simplerockethody  is  modeled  with  50  flat  facets.  The  curved  truncated  conic 
sections  of  lUSrocketbody  and  SL12rocketbody  are  modeled  with  20  flat  facets.  To  produce  accurate  light  curves  of 
curved  surfaces,  however,  requires  a  larger  number  of  facets  depending  on  the  BRDF  exponential  factor  of  the  sur¬ 
face.  This  is  discussed  further  in  the  appendix. 


4.  SOLAR  RADIATION  PRESSURE  (SRP)  MODEL 

The  acceleration  caused  by  the  SRP  can  be  calculated  by  summing  the  individual  contributions  of  all  the  constit¬ 
uent  illuminated  “facets”  that  make  up  the  object.  The  acceleration  is 
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where  Fi{X)  is  the  solar  flux  (in  Wm'^nm'^),  Ak  is  the  facet  area,  fu  is  the  fraction  of  the  facet  that  is  illuminated  (due 
to  self- shadowing),  niso  is  the  mass  of  the  object,  c  is  the  speed  of  light,  and  the  BRDF  for  each  facet  is  integrated 
over  all  observer  directions  and  all  wavelengths.  Additionally,  (v)+  =  xH{x)  where  H{x)  is  the  Heaviside  step  func¬ 
tion  which  is  unity  for  positive  values  and  zero  elsewhere. 
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For  certain  BRDFs,  the  integral  can  be  solved  analytically.  For  example,  the  BRDF  with  a  Lambertian  diffuse 
component  and  purely  “mirror-like”  specular  component  (where  R  =  2(l-n)n-L  is  the  direction  of  mirror-like 

reflection  of  L )  is  given  by 


fr=d 


+  s\ 


v(v-r)~ 
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(8) 


where  6  is  the  Delta  function,  yields  an  acceleration  due  to  the  SRP  of 
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where  Fsun  is  the  total  solar  flux  over  all  wavelengths.  This  is  the  conventional  formula  typically  used  to  calculate 
the  SRP  acceleration  for  a  facet-based  model  (e.g.  equivalent  to  equations  derived  in  Appendix  A  in  Ref.  [5]). 

For  a  more  complicated  BRDF,  such  as  the  Ashikhmin- Shirley  model  described  previously,  the  exact  solution 
obtained  by  numerically  integrating  Eq.  (7)  is  different  than  the  idealized  solution  obtained  by  Eq.  (9).  Ref.  [1]  de¬ 
veloped  a  way  of  computing  correction  factors  for  Eq.  (9)  and  the  acceleration  is  calculated  by  using 
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where  the  A  factors  depend  on  the  illumination  angle  and  parameters  within  the  BRDF  model.  Thus,  for  Aj  =  1,  the 
diffuse  reflectance  is  Lambertian  and  for  A^i  =  A^2  =  1,  the  specular  reflectance  is  mirror-like.  The  further  these  A 
factors  are  from  unity,  the  greater  the  difference  between  the  BRDF  model  and  the  idealized  BRDF  of  Eq.  (9). 

The  SRP  acceleration  is  highly  dependent  on  the  object’s  attitude.  For  the  rocket  body  models,  it  is  assumed 
they  are  single  axis  rotators  about  the  model’s  x-axis  and  have  one  of  three  initial  attitudes  (corresponding  to  the 
rotation  axis  in  the  x,  y,  and  z  ECI  directions)  and  one  of  two  initial  angular  rates  (corresponding  to  a  “slow”  rota¬ 
tion  with  a  period  of  about  1  hour  and  a  “fast”  rotation  with  a  period  of  about  5  minutes).  For  the  debris  models,  the 
same  initial  attitude  is  chosen  for  all  cases  (corresponding  to  body  frame  equal  to  ECI  frame)  and  one  of  four  initial 
angular  rates  (corresponding  to  two  different  “slow”  tumbles  and  two  different  “fast”  tumbles  where  “slow”  and 
“fast”  are  comparable  to  the  rocket  body  angular  rates).  For  all  cases,  the  SRP  model  of  Eq.  (9)  and  Eq.  (10)  are 
also  simulated. 


5.  SIMULATIONS  SETUP  AND  RESULTS 

The  orbit  was  initialized  to  that  of  the  two-line  element  (TLE)  of  the  near  geosynchronous  rocket  body  SSN 
19550  with  observations  simulated  starting  at  2013  May  01  at  0:00:00  UT  and  extending  to  2013  Aug  01  at  0:00:00 
UT  (92  days).  Specifically,  the  initial  Cartesian  state  is  r  =  [-41557.77598  -6572.93284  1836.98724]  km  and  v  = 
[0.43217  -2.94731  -0.76895]  km/s  corresponding  to  an  orbit  with  an  eccentricity  of  0.001568  and  inclination  of 
14.6755  degrees. 

For  the  rocket  body  models,  there  are  a  total  of  144  simulated  cases  for  each  of  the  three  shape  models.  For  the 
debris  models,  there  are  a  total  of  96  simulated  cases  for  each  of  the  three  shape  models.  Figures  3  and  4  display 
the  total,  radial,  in-track,  and  cross-track  distances  from  the  baseline  of  no  SRP  acceleration  for  all  cases  corre¬ 
sponding  to  simplerocketbody  and  simpledebrismodel(  1 )  respectively. 
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Table  2.  Simulation  Results 


name 

AMR 

(max) 

Average  total 
distance  (km) 

Range  in  total 
distance  (km) 

simplerocketbody 

0.0180 

28.4 

11.6 

lUSrocketbody 

0.0115 

15.0 

8.3 

SLl  2  rocketbody 

0.0076 

8.1 

2.6 

simpledebris(  1 ) 

2.0 

1215.3 

1695.4 

simpledebris(2) 

0.4 

145.8 

322.4 

crinkleddebris 

0.4 

131.3 

369.8 

Table  2  lists  the  average  total  distance  and  spread  in  total  distances  of  the  full  set  of  cases  from  the  baseline  after 
the  92  day  propagation.  In  both  the  rocket  body  and  debris  simulations,  there  is  a  clear  direct  relationship  between 
the  maximum  area-to-mass  ratio  (AMR)  and  both  the  average  total  distance  from  baseline  and  the  range  in  total  dis¬ 
tances  from  baseline  after  92  days.  This  is  shown  in  Fig.  5. 
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Figure  5.  Average  total  distance  and  spread  in  total  distances  from  baseline  after  92-day 
propagation  for  the  six  shape  models 

The  total  difference  from  a  baseline  while  varying  a  particular  surface  or  environmental  parameter  can  be  calcu¬ 
lated  using  all  cases  corresponding  to  each  of  the  six  shape  models.  For  example,  using  the  lowest  value  for  the 
BRDF  parameter  as  the  baseline.  Figs.  6  and  7  show  the  difference  of  the  p  =  Fo  =  0.5  (red)  cases  and  p  =  Fo  =  0.9 
(green)  cases  from  the  p  =  Fo  =  O.l  cases  where  all  the  other  parameters  are  the  same;  Figs.  8  and  9  show  the  differ¬ 
ence  of  the  J  =  0.8  cases  from  the  d  =  0.2  cases  where  all  the  other  parameters  are  the  same;  and  Figs.  10  and  11 
show  the  difference  of  the  ^  =  500  cases  from  the  ^  =  50  cases  where  all  the  other  parameters  are  the  same. 
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Figure  6.  Orbital  differences  of  rocket  body  models  for  BRDF  reflectance  parameter 


Figure  7.  Orbital  differences  of  debris  models  for  BRDF  reflectance  parameter 


Figure  8.  Orbital  differences  of  rocket  body  models  for  BRDF  diffuse  fraction  parameter 
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simpledebris(l)  simpledebris(2)  crinkleddebris 
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Figure  9.  Orbital  differences  of  debris  models  for  BRDF  diffuse  fraction  parameter 

simplerocketbody  lUSrocketbody  SL12rocketbody 


time  (days)  time  (days)  time  (days) 

Figure  10.  Orbital  differences  of  rocket  body  models  for  BRDF  exponential  factor  parame¬ 
ter 

simpledebris(l)  simpledebris(2)  crinkleddebris 
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Figure  11.  Orbital  differences  of  debris  models  for  BRDF  exponential  factor  parameter 

Figs.  12  and  13  display  the  total  difference  of  the  simplified  SRP  model  from  the  physically-consistent  SRP 
model  with  all  other  parameters  remaining  the  same. 
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Figure  12.  Orbital  differences  of  rocket  body  models  for  SRP  model  used 


Figure  13.  Orbital  differences  of  debris  models  for  SRP  model  used 

The  same  can  be  done  for  the  attitude  and  angular  rates  related  cases.  The  difference  of  the  initial  rotation  axis 
along  the  y-axis  ECI  cases  and  initial  rotation  axis  along  the  z-axis  ECI  cases  from  the  initial  rotation  axis  along  the 
x-axis  ECI  cases  where  all  the  other  parameters  are  the  same;  the  difference  of  the  “fast”  initial  angular  rate  cases 
from  the  “slow”  initial  angular  rate  cases  where  all  the  other  parameters  are  the  same;  and  the  difference  of  the  two 
“fast”  initial  angular  rate  cases  and  the  “slow”  initial  angular  rate  cases  from  the  other  “slow”  initial  angular  rate 
cases  where  all  the  other  parameters  are  the  same.  Tables  3  and  4  summarize  the  results  by  tabulating  the  average 
distance  and  maximum  distance  from  baseline  for  each  of  the  changing  parameters  and  for  each  of  the  rocket  body 
and  debris  models. 


Table  3.  Average  Distance  (km)  from  Baseline 


name 

P 

d 

n 

Initial 

attitude 

Initial  rate 

SRP  model 

simplerocketbody 

2.98 

3.39 

0.22 

3.62 

0.23 

0.46 

lUSrocketbody 

1.99 

1.72 

0.16 

2.36 

0.72 

0.36 

SL12rocketbody 

0.82 

0.88 

0.07 

0.37 

0.08 

0.16 

simpledebris(  1 ) 

198.1 

284.1 

63.0 

- 

561.6 

151.2 

simpledebris(2) 

45.7 

64.6 

18.5 

- 

109.3 

34.3 

crinkleddebris 

70.9 

78.5 

40.0 

- 

115.8 

71.4 
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Table  4.  Maximum  Distance  (km)  from  Baseline 


name 

P 

d 

n 

Initial 

attitude 

Initial  rate 

SRP 

model 

simplerocketbody 

7.25 

9.35 

1.12 

6.13 

2.15 

1.63 

lUSrocketbody 

5.03 

4.95 

1.59 

5.83 

1.90 

2.01 

SL12rocketbody 

1.92 

2.15 

0.35 

0.90 

0.36 

0.41 

simpledebris(  1 ) 

880.9 

1359.0 

805.0 

- 

1346.0 

613.3 

simpledebris(2) 

297.8 

289.3 

310.5 

- 

298.5 

320.5 

crinkleddebris 

319.4 

417.3 

357.1 

- 

456.0 

321.5 

6.  CONCLUSION 

Whereas  the  choice  of  solar  radiation  pressure  (SRP)  model  (whether  physically  consistent  or  simplified)  does  have 
a  measureable  effect  on  the  perturbations  of  rocket  body  and  debris  objects  in  near-geosynchronous  orbits,  it  is  a  less 
significant  than  changes  in  other  parameters  that  govern  the  magnitude  of  the  SRP  acceleration.  Namely,  these  pa¬ 
rameters  include  those  associated  with  the  bidirectional  reflectance  distribution  function  (BRDF)  (reflectance,  dif¬ 
fuse  fraction,  and  exponential  factor)  and  attitude  (initial  attitude  and  angular  rates).  Most  important  of  these  appear 
to  be  the  overall  diffuse  and  specular  reflectance  and  the  diffuse  to  specular  fraction.  In  the  rocket  body  cases  where 
the  object  started  as  a  simple  rotator,  the  initial  attitude  was  also  relatively  significant  while  the  exponential  factor 
and  initial  angular  rate  along  with  the  SRP  model  were  less  so.  In  the  debris  cases  where  the  object  started  in  a  tum¬ 
bling  state  (and  initial  attitude  was  not  changed),  the  initial  angular  rates  appeared  to  also  be  relatively  significant 
while  only  the  exponential  factor  and  SRP  model  were  less  so. 
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